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A theoretical analysis of the anomalous diffusion transport mechanism suggests a possible con- 
nection between wall current drain and magnetic flux through the orbital trajectories of charged 
particles in a plasma submitted to a strong magnetic field. Then Bohm diffusion coefficient is 
straightforwardly obtained. 

PACS numbers: 05.70.-a,52.25.Fi,52.55.-s,55.55.Dy,52.40.Hf 



The problem of the plasma-wall interactions is of fun- 
damental importance in plasma physics. Studying the 
anomalously high diffusion of ions across magnetic field 
lines in Calutron ion sources (electromagnetic separator 
used by E. Lawrence for uranium isotopes) gave the firsts 
indications of the onset of a new mechanism bimon 
advanced first with a suggestion that the observed losses 
could be explained by the highly anisotropic medium in- 
duced by the magnetic field lines, favoring electron cur- 
rent to the conducting walls - the electron "short-circuit" 
problem Q. Experiments done by Geissler Q in the 
1960's have shown that diffusion in a plasma across a 
magnetic field was nearly classical (standard) diffusion 
when insulating walls impose plasma ambipolarity, but 
in the presence of conducting walls charged particles dif- 
fuse at a much higher rate. 

This problem of plasma-wall interaction becomes more 
complex when a complete description is aimed of a mag- 
netized nonisothermal plasma transport in a conducting 
vessel 0- In the area of fusion reactors, there is strong 
indication that for plasmas large but finite Bohm-like dif- 
fusion coefficient appears above a certain range of B 

Progress in the understanding of the generation of con- 
finement states in a plasma is fundamental |6j to pur- 
sue the dream of a fusion reactor 0, 0- Anomalous 
diffusion is a cornerstone in this quest, as recent re- 
search with tokamaks suggest that the containment time 
t w 10 8 i? 2 /2Db, with R denoting the minor radius of 
a tokamak plasma and Db is the Bohm diffusion coef- 
ficient 0- Controlled nuclear fusion experiments have 
shown that transport of energy and particles across mag- 
netic field lines is anomalously large (i.e, not predicted 
by classical collision theory). 

The conjecture made by Bohm is that the diffusion 
coefficient is Db = akT/eB, where T is the plasma 
temperature and a is a numerical coefficient, empiri- 
cally taken to be 1/16 Q. Usually the origin of the 
anomalous diffusion has been assumed to be due to the 
turbulence of small-scale instabilities (see, for example, 
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Refs. HOHEl)- 

In a seminal paper |12| a conjecture was proposed 
based on the principle of minimum entropy-production 
rate, stating that a plasma will be more stable whenever 
the internal product of the current density by an elemen- 
tary conducting area dA at every point of the boundary 
- excluding the surface collecting the driving current - is 
null, (j • A) = 0. 

The work reported in this Letter proposes a mech- 
anism of wall current drain set up together with the 
frozen-in effect of the magnetic field as a possible alterna- 
tive explanation of the anomalous diffusion mechanism. 
From collisional low temperature plasmas to a burning 
fusion plasma subject the plasma confinement vessel to 
strong wall load, both in stellarator or tokamak operating 
modes, this explanation could be of considerable interest. 

The general idea proposed by Robertson ^3] assume 
that the plasma boundary is composed of small elements 
of area Ai, each one isolated from the others but each one 
connected to the exterior circuit through its own resistor 
Ri and voltage Vi. The entropy production rate in the 
external circuits is 

§ = £^-^) 2i? " (i) 

i 

where T is the temperature of the resistors, supposed to 
be in thermal equilibrium with all the others. 

In fact, a straightforward application of Eq.^to a cold 
plasma made of electrons and just one ion component 
gives 

dS e 2 

— = — {~n e u e E+D e Vn e + n l u i I]-D l Vn i ) 2 A 2 R. (2) 
at 1 

Under the usual assumptions of quasi-neutrality and 
quasi-stationary plasma (see, for example, Ref. p^|) 

^ = e = const.; n e v e = rijVj. (3) 
Hence, Eq. 2] becomes 
dS e 2 

— = — [E(e/ij - u e )n e + Wn e (D e - D l e)] 2 A 2 R. (4) 
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If there is no entropy production S = 
ambipolar electric field is recovered ^ 



E 



[lit n e 



0, and then the 



(5) 



This conceptual formulation provides new insight regard- 
ing ambipolar diffusion. In a thermal equilibrium state, a 
plasma confined by insulating walls will have an effective 
coefficient given by the above Eq.JS] a situation frequently 
encountered in industrial applications. This example by 
itself relates ambipolar diffusion with no entropy pro- 
duction in the plasma. Allowing plasma currents to the 
walls, entropy production is greatly enhanced, generating 
altogether instabilities and plasma losses ^1 • 

But confined plasmas are in a far-nonequilibrium state 
(with external surroundings) and it is necessary to es- 
tablish a generalized principle that rule matter, which 
we develop in the next lines. 

Experiments give evidence of transport of particles and 
energy to the walls At the end of the 1960s, ex- 

perimental results obtained in weakly ionized plasma |3( 
and in a hot electron plasma |l5j (this one proposing 
a possible mechanism of flute instability) indicated the 
strong influence conducting walls have on plasma losses 
across magnetic field lines. Geissler suggested that 
the most probable explanation was due to the existence 
of diffusion-driven current flow through the plasma to the 
walls. Beilinson et al. Q shown the possibility to control 
the profiles of plasma parameters by applying potential 
difference to various parts of the conducting walls. Con- 
cerning fusion reactions, Taylor [lOj provided a new in- 
terpretation of tokamak fluctuations as due to an inward 
particle flux resulting from the onset of filamentary cur- 
rents. 

We consider a simple axisymmetric magnetic configu- 
ration with magnetic field lines parallel to z-axis with a 
plasma confined between two electrodes (see Fig. I). In 
general terms, a particle motion in a plasma results in 
a massive flux. As long as the flux is installed, the flux 
will depends naturally on a force F - in this case the 
pressure gradient-driven process of diffusion to the wall 
- responsible of the wall driven current j. According to 
the fundamental thermodynamics relation, the plasma 
internal energy variation dU is related to the amount of 
entropy supplied or rejected and the work done by the 
driven force, through the equation 



dU 
~dt 



= (j • A) R + F • 



dr 
~dl 



(6) 



The last term we identify with the macroscopic diffusion 
velocity depicting the process of plasma expansion to 
the wall. To simplify somehow the calculations we as- 
sume a single plasma fluid under the action of a pressure 
gradient (F = ALdp/dy). 

In the presence of steady and uniform magnetic field 
lines (this simplifies the equations, but do not limit the 
applicability of the model), the particles stream freely 



along them. From magnetohydrodynamic we have a kind 
of generalized Ohm's law (see, for example, Ref. 16]) 



Vp = -enE - en[v x B] + [j x B] 



en] 
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(7) 



where a = e 2 nT e /m e is the electric conductivity, with 
r e denoting the average collision time between electrons 
and ions. Under the equilibrium condition 



Vp= [j x B], 



(8) 



holds. Therefore, after inserting Eq. [S] into Eq. [7| it is 
obtained the y component of velocity 



E x 
~B 



1 dp 



(9) 



From Eq.OHwe have the classical diffusion coefficient scal- 
ing with I / B 2 and thus implying a random walk of step 
length Th (Larmor radius). To get the anomalous diffu- 
sion coefficient we must consider the process of diffusion 
to the wall - in the presence of an entropy source - with 
the combined action of the wall current drain, as already 
introduced in Eq. 

Therefore, using the guiding center plasma model the 
particle motion is made with velocity given by 



enw d 



[Vp x B] 
B 2 ■ 



(10) 



This equation form the base of a simplified theory of 
magnetic confinement. In fact, the validity of Eq. ITUI 
is restrained to the high magnetic field limit, when the 
Larmor radius is shorter than the Debye radius. 

Considering motion along only one direction perpen- 
dicular to the wall (y-axis), it is clear that 



A 2 [dp 



(id 



If we consider a quasi-steady state plasma operation, 
the plasma total energy should be sustained. Hence, 
dU j dt = 0, and the power associated with the driven 
pressure-gradient is just maintaining the dissipative pro- 
cess of plasma losses on the wall. Eq. [H] govern the evo- 
lution of diffusion velocity. Hence, we have 



nv d = 



nRA kT dn 
L TPdy 



„ dn 
-D T —, 
ay 



(12) 



with Dt denoting the transverse (across the magnetic 
field) diffusion coefficient given by 



D T = 



nRA kT 
L B 2 ' 



(13) 



This new result coincides with the classical diffusion coef- 
ficient |l3| | whenever nRA/ L = mis e i/e 2 , containing a de- 
pendence on collision frequency and particle number den- 
sity. Others theoretical approaches to this problem were 
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advanced by Bohm |l|, who proposed an empirically- 
driven diffusion coefficient associating plasma oscillations 
as the source of the enhanced diffusion, while Tonks 
have shown that the current density that is present in 
a magnetically immobilized plasma is only generated by 
the particle density gradient, not being associated with 
any drift of matter. Simon electron "short-circuit" 
proposes an explanation for the different rates of diffu- 
sion electrons and ions do experiment across the mag- 
netic field as due to an unbalance of currents flowing to 
the wall. 

In the absence of collisions, the guiding centers of 
charged particles behave as permanently attached to the 
same lines of force. On the contrary, as a result of col- 
lisions with others charged particles the guiding centers 
shift from one line of force to another resulting in a dif- 
fusion of plasma across the field lines. In our model, 
each orbit constitutes an elementary current / eventu- 
ally crossing the wall. 

However, the particle diffusion coefficient as shown in 
Eq. 1131 gives evidence of an interplay between the resis- 
tance the elementary circuit is submitted when in contact 
with the walls in the presence of the frozen-in effect. In 
fact, for sufficiently strong magnetic fields apparently a 
hydrodynamic behavior of the plasma is installed [filllSj. 
with the appearance of "convective cells" and the 1/B 
behavior dominates, giving birth to the anomalous dif- 
fusion mechanism. The onset of freezing magnetic lines 
is valid whenever the Lundquist number S ^> 1 (convec- 
tion of the magnetic field dominated medium). In this 
case the magnetic field lines are frozen-in in the medium 
(consequence of a vortex type of character of the mag- 
netic field B) and the flux of them across a given surface 
is constant: 

$ = BA' = BL 2 a. (14) 

Remark that A' is now the surface delimited by the ele- 
mentary circuit 7 (see Fig.^) and a < 1 is just a geomet- 
rical factor (e.g. a — 7r/4 at the limit of a circular orbit). 
This situation is fundamental to the onset of anomalous 
diffusion. Free electrons orbits are helical, but as Fig. 
shows, their projections at right angles to the field are 
circular. Each particle orbit constitute an elementary 
circuit with B-field cutting its surface being associated 
with it an elementary flux $. At the same time we can 
envisage each orbit as constituting by itself an elementary 
circuit, some of them intersecting the wall and thus the 
circuit is closed inside the wall. Therefore a resistance 
R drags the charged flow at the conducting wall. It is 
therefore plausible to associate to this elementary circuit 
a potential drop V and the all process being equivalent 
to a current / flowing through the elementary circuit. 

Assuming the plasma is a typical weakly coupled, hot 
diffuse plasma with a plasma parameter (number of par- 
ticles in Debye sphere) A = n\ 5 De 1, it is more likely to 
expect nearly equal average kinetic and potential energy. 
However, the typical plasma parameter encountered in 
glow discharges or in nuclear fusion is A 3> 1. This means 
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FIG. 1: Schematic of the geometry for the plasma-wall cur- 
rent drain model. The uniform magnetic field points down- 
ward along Oz. Particles describe orbits in the plane xOy 
intersecting the wall (plan xOz). Orbits are represented by a 
semi-circular line for convenience. L is the maximum distance 
the trajectory attains from the wall. 

that the average kinetic energy is larger than the average 
potential energy. To contemplate all range of A we can 
relate them through the relationship 

pV = (J ■ A)S. (15) 

Here, p is the charge density, A is the vector potential, 
J is the current density and 5 < 1 is just a parameter 
representing the ratio of potential to kinetic energy. Of 
course, when A > 1, then 5 < 1. This basic assumption is 
consistent with the hydrodynamic approximation taken 
in the development of equations. The limitations of the 
model are related with the unknowns A and S that can 
be uncovered only through a self-consistent model of the 
plasma. However, our analysis of anomalous diffusion re- 
mains general and added new insight to the phenomena. 

Now suppose that the diffusion current is along y-axis 
J = —Jyiiy (see Fig.l). Consequently, A = — A y u y , and 
then the potential drop will depend on x-coordinate: 

P [V{ Xl ) - V(x )} = J y [A y { Xl ) - A y (x )]S. (16) 

Multiplying both members by the area A' — x\Z\ and 
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length L = yi, we have 

QAV = Iyi[A y (xx) - A y (x )]S = I<PS. (17) 

$ = / 7 (A • g?x) is the flux of the magnetic field through 
the closed surface bounded by the line element cfoc (ele- 
mentary circuit 7, see also Fig^). By other side, natu- 
rally, the total charge present on the volume V = XiyiZi 
is such as Q — ie, with i an integer. This integer must be 
related to ions charge number. From Eg. 1171 we obtain 



R 



AV 



r $ ,BL 2 
o— = ao— — . 



(18) 



But, the particle density is given by n = N/LA, with N 
being now the total number of charged particles present 
in volume V = AL. Since i — N, we retrieve finally the 
so-called Bohm-diffusion coefficient 



D, 



aS- 



.kT 
) — . 

eB 



(19) 



However, this expression suffers of the indetermination 
of the geometrical factor a. This factor is related to 
the ions charge number, it depends on the magnetic field 
magnitude and as well the external operating conditions 
(due to increased collisional processes, for ex.). The exact 
value of the product ad can only be determined through a 
self-consistent plasma model, but we should expect from 
the above discussion that aS < 1. Furthermore, Eg. 1181 
can be used as a boundary condition (simulating an elec- 
trically floating surface) imposed when solving Poisson 
equation. 

Also it worth to emphasize that when inserting Eq. IT51 
into Eq. 1131 and considering the usual definition of mo- 
mentum transfer cross section, then it can be obtained 
a new expression for the classical diffusion coefficient as 
a function of the ratio of collisional v and cyclotron fre- 
quency f2, although (and in contrast with the standard 
expression), now also dependent on the geometrical fac- 
tor a and energy ratio 5: 



Dn 



(aS) 



v kT 
f2 m 



(20) 



This explains the strong dependence of the classical dif- 
fusion coefficient on is/tt showing signs of anomalous dif- 
fusion as discussed in Ref. [HJ (obtained with a time 
resolved Monte Carlo simulation in an infinite gas under 
uniform fields) and, in addition, the strong oscillations 
shown up in the calculations of the time dependence of 
the transverse component of the diffusion tensor for elec- 
trons in low-temperature rf argon plasma. Those ba- 
sic features result on one side from its dependence on 
R, which is proportional to the flux. Therefore, a flux 
variation can give an equivalent effect to the previously 
proposed mechanism: whenever a decrease (or increase) 
in the flux is onset through time dependence of electric 
and magnetic fields, it occurs a strong increase (or de- 
crease) of the diffusion coefficient. By other side, when 
the resistance increases it occurs a related decrease of 
charged particles tangential velocity and its mean energy. 
So far, this model gives a new insight into the results re- 
ferred in ^3 and also it explains why the same effect is 
not obtained from the solution of the non-conservative 
Boltzmann equation as applied to an oxygen magnetron 
discharges with constant electric and magnetic fields [2(| . 



To summarize, we introduced in this Letter a simple 
mechanism providing an interpretation of the anomalous 
diffusion in a magnetized confined plasma. In fact, above 
a certain range of B the magnetic field frozen-in effect 
is settled in the plasma, implying a constant magnetic 
field flux through the elementary orbits of the charge 
carriers. Whenever conducting walls are bounding the 
plasma current drain to the walls occurs naturally and a 
Bohm-like behavior of the transverse diffusion coefficient 
results. The suggested mechanism could lead to a better 
understanding of the mechanism of plasma-wall interac- 
tion and help to develop a full-scale numerical modelling 
of present fusion devices or collisional low-temperature 
plasmas. 
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